Cuprate Fermi orbits and Fermi arcs: the effect of short-range antiferromagnetic order 
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We consider the effect of a short antiferromagnetic correlation length £ on the electronic band- 
structure of the underdoped cuprates. Starting with a Fermi-surface topology consistent with that 
detected in magnetic-quantum-oscillation experiments, we show that a reduced £ gives an asym- 
metric broadening of the quasiparticle dispersion, resulting in simulated ARPES data very similar 
to those observed in experiment. Predicted features include the presence of 'Fermi arcs' close to 
ak = (7r/2,7r/2), where a is the in-plane lattice parameter, without the need to invoke a d-wave 
pseudogap order parameter. The statistical variation in the fc-space areas of the reconstructed 
Fermi-surface pockets causes the quantum oscillations to be strongly damped, even in very strong 
magnetic fields, in agreement with experiment. 



Magnetic quantum oscillations such as the de Haas- 
van Alphen and Shubnikov-de Haas effects are acknowl- 
edged to be the most reliable method to determine the de- 
tailed low-temperature Fermi-surface topologies of met- 
als P, 0] • Recently, Shubnikov-de Haas oscillations have 
been observed in the underdoped cuprate superconduc- 
tors YBa 2 Cu 3 6 .5 [3] and YBa 2 Cu 4 8 @, Hf; the data 
reveal quasi-two-dimensional Fermi-surface sections with 
cross-sections that are ~ 2 — 2.4% of the area of the 
Brillouin zone. This result is consistent with a pic- 
ture [3J] in which antiferromagnetic order breaks transla- 
tional symmetry, reconstructing the large Fermi surface 
seen in overdoped cuprates |6j into four small, symmetry- 
related pockets. At first sight, this conclusion is at vari- 
ance with angle-resolved photo-emission (ARPES) ex- 
periments 0, Q , which find only so-called "Fermi arcs" . 
However, we show here that the quantum-oscillation and 
ARPES results are compatible. Starting with a Fermi 
surface consistent with the quantum-oscillation experi- 
ments, we show that a reduced antiferromagnetic corre- 
lation length £ gives an asymmetric broadening of the 
quasiparticle dispersion, resulting in predicted ARPES 
spectral weights similar to those observed in experiment. 

The breaking of translational symmetry due to anti- 
ferromagnetism can be described by a modulation vec- 
tor Q that is sharply defined in an ideal antiferromag- 
net [9(; this is the case in the parent Mott-insulating 
phase of the cuprates [1(3], realized for hole dopings 
p < 0.05. However, in the underdoped regime (0.05 < 
p < 0.2), slow antiferromagnetic fluctuations, detected 
in neutron-scattering experiments, cause Q to become 
subject to statistical variations over a wide range of con- 
ditions [HI EH]- At temperatures T < 10 K, the fluc- 
tuations become sufficiently slow for the system to re- 
semble a glassy state comprising antiferromagnetic mi- 



crodomains with lifetimes Ti ~ 10 s [12]]. This life- 
time exceeds by several orders of magnitude both that 
for the re-equilibration of the Fermi sea occurring in 
ARPES experiments (r r ~ 10" 15 s) @, @] and the pe- 
riod of quasiparticle Fermi-surface orbits in strong mag- 




netic fields (27rw c 



10" 



s, where u c is the cyclotron 



frequency 0, H|)- Therefore, from the perspective of 
both experimental techniques, quasiparticles can be con- 
sidered to propagate through a glassy antiferromagnetic 
medium in which the staggered magnetization is qua- 
sistatic and periodic over distances of order £. 
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FIG. 1: Quadrant of a notional cuprate Fermi surface for 
p — 0.095 doped holes per planar Cu atom. a. The Fermi 
surface given by Eq. [1] in the absence of ordering, b. The 
Fermi surface (white line) after reconstruction (i.e. folding 
about the dotted line) with respect to aQ = aQo = (n, n) 
with A = tio in Eq. [2] The red line shows the corresponding 
Fermi surface for a slightly different Q, shifted relative to Qo. 

To see how ARPES and quantum oscillation experi- 
ments are affected by the short antiferromagnetic correla- 
tion lengths £, we consider a simplified, two-dimensional 
quasiparticle dispersion of the form 

£k = 2tio(cosafca; + cos ak y ) + 

2tn(cosa(k x + k y ) + cosa(k x — k y )). (1) 

Eq. Q] gives rise to a large hole Fermi surface sheet cen- 
tered on ak = (n,Tr) for t 10 < 0, tu ~ — 0.34x< 10 and 
p > 0; this is approximately representative of cuprate 
Fermi surfaces in the absence of ordering [|| . Figure [T^, 
shows a quadrant of such a Fermi surface for p — 0.095 
holes per Cu atom within the Cu-0 planes. The entire 
Fermi surface accommodates 1 + p holes per unit cell (in 
the absence of ordering), having a fc-space cross-sectional 
area = 27r 2 (l +p)/a 2 , where a is the in-plane lattice 
constant. The Fermi energy ep is obtained self consis- 
tently from Ak — J J f(£F — ek)dk x dk y , where f(x) is the 
Fermi-Dirac distribution function with the double inte- 
gration being performed within the first Brillouin zone. 
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The reconstruction of such a Fermi surface caused by 
long-range antiferromagnetic order is well understood [9j . 
The staggered magnetization m = m cos(r • Q) repre- 
sents a new periodicity of the system that introduces a 
new, smaller Brillouin zone and causes the quasiparticle 
dispersion relationship (Eq. [Q) to become gapped close 
to the zone boundaries by an energy 2 A. In the weak- 
coupling limit, A <C 4£io, the modified dispersion can be 
written 



e k,Q 



; k+Q 



± 



(2) 



in the absence of an applied magnetic field [l3| ■ The 
white line in Fig. [TJd shows one of the four equivalent 
Fermi-surface pockets that result for A = t\o and a|Q| = 
a IQo| = |(tT) ti") I To accommodate p holes within the 
reduced Brillouin zone (bounded by the dotted line), each 



pocket has an area Ak,Q = n 2 p/a 2 [15]. To conserve 
charge neutrality, the introduction of antiferromagnetic 
order is accompanied by a shift in Fermi energy from £p 
to £f,q ; hence A k , Qo = \ j j f(e F ,Q - £k,Q )dfc a; dfc y . 

The parameters p — 0.095 and A = tio used in Fig. [TJd 
are chosen to mimic the experimental Fermi surfaces of 
Cai. 9 Nao.iCu0 2 Cl2 and YBa 2 Cu 4 8 [4] for which 
the corresponding quantum-oscillation frequency is F — 
hp/Aea 2 « 650 T. Whereas the Mott insulating parent 
phase (p < 0.05) may be understood within a strong- 
coupling description of antiferromagnetism 2, 17 1, the 
appropriate regime for p w 0.1 remains a matter of con- 



jecture [11J. In the strong coupling limit (A ^> t w ), the 



hole pockets at (ir/2,ir/2) would depart only marginally 
from the form depicted in Fig. [T|d, the primary effect 
being to renormalize the widths of the upper and lower 
bands defined by Eq. [2] In the very weak coupling limit 
(A/tio — > 0), however, in addition to the hole pockets 
at (tt/2, 7t/2), small electron pockets would appear at 
(■7T, 0); there is as yet no evidence for the latter pockets 
in quantum-oscillation data jij. 

We model the effect of short correlation lengths using 
a probability distribution for Q |18], as determined by 
neutron-scattering experiments [lfj. As stated above, in 
ARPES and quantum-oscillation experiments we are in 
the limit re ^> a;" 1 3> r r , so that the distribution can be 



considered static: 



K 



PQ 



r 2 + (Qo - q) 



(3) 



with J J pQdQxdQy = 1 [13]. To help understand how 
the Fermi surface is modified by variations in Q, the red 
curve in Fig. [TJd shows how the pocket is deformed by a 
small shift of Q relative to Qo- Whereas the 'A' branch 
of the Fermi-surface orbit (required by the translational 
symmetry of the new Brillouin zones) is strongly modi- 
fied by this shift, the 'B' branch that corresponds more 
closely to the original Fermi surface is less altered. This is 



true irrespective of the direction of the shift in Q. Thus, 
the 'B' branch is more resilient against statistical varia- 
tions in Q than is 'A'. On considering all possible values 
of Q within this distribution, the resulting dispersion and 
Fermi surface become broadened; observables such as the 
photoemission spectral weight and Landau-level broad- 
ening can then be estimated by analogy with other types 
of "phase smearing" [3, [5(| • 

ARPES provides information on £k and the extent to 
which it is broadened by correlations or other factors, giv- 
ing rise to a photoelectron intensity that is distributed in 
energy and fc-space @, H, Hl[. When the energy distri- 
bution curve is cut off by the Fermi-Dirac distribution, 
the rapidity with which the integrated intensity grows 
with Ep — £k close to the Fermi energy can be used to 
obtain a 'Fermi surface intensity plot' [3, H, EH- The dis- 
continuous jump at £f encountered in conventional met- 
als therefore gives rise to an associated spectral weight 
Sk = ^( £ k — £f) that is infinite at T = under ideal con- 
ditions. In actual experiments, however, with overdoped 
Tl2Ba2Cu06+a providing a good example, the spectral 
weight at £p is broadened by finite temperature and self- 
energy effects and by instrumental energy resolution 21 1 . 

In the case of antiferromagnetism with a distribution of 
Q, we determine the spectral weight from the convolution 
of <5(£k,Q - £f,q ) with Eq. [3j 



Sk.Qo = J J <5(£k,Q - £F,Q )PQdQ x dQy 



(4) 



The sensitivity of the reconstructed Fermi surface to 
statistical variations in Q (i.e. Fig. [Qa) introduces a 
strongly fc-space dependent broadening of the spectral 
weight. Fig. [2] shows the spectral weight calculated us- 
ing Eq. [f] for several different values of £. At large £, 
the spectral weight strongly resembles the outline of the 
Fermi surface depicted in white in Fig.[T)3. As £ becomes 
shorter, however, the spectral weight corresponding to 
the part of the Fermi surface produced by the fc-space 
reconstruction (labeled 'A') fades to be replaced by a 
greater concentration along an arc of the original unre- 
constructed Fermi surface. At £ = 10 and 20 A, the 
distribution of spectral weight closely resembles that ob- 
served in Cai.9Nao.iCu02Cl2 (see Figs. lb,c of Ref. 0]), 
for which p w 0.1. These values of £ are in excellent 
agreement with neutron-scattering estimates made at the 
same nominal p in a variety of cuprates (at temperatures 
~ tens of kelvin in zero applied magnetic field) [Toj . 

The model reproduces the features seen in 
Cai. 9 Nao.iCu0 2 Cl2 and Bi 2 Sr 2 CaCu20 8+< 5 ARPES 
data @, H|]. Fig. [3h. shows the angle (9) dependence 
of the spectral weight from Fig. [2] (£ = 10 A) tracing 
the k = kp arc of the original unreconstructed Fermi 
surface in Fig. QJi. The sinusoidal variation of Sk F ,Q 
with 9 (Fig. [3k,) and the strongly reduced gradient of 
the energy distribution intensity / with respect to £ at 
k ±45 o in Fig. [3b are very similar to ARPES data @, [!]• 
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FIG. 2: Spectral weight calculated over a quadrant of the 
Fermi surface using Eq. [4] for different values of £. Since the 
continuous variables k and e are replaced by discrete quan- 
tities in the simulations, the delta function is replaced by a 
Kronecker delta. 



Such features have been interpreted in terms of a 
'pseudogap'; note that our model shows that short-range 
antifcrromagnetic order can also give the 'illusion' of 
a rf-wave order parameter with a node at 6 = 0° and 
antinodes at 8 = ±45° (see inset) Q- This is also shown 
in Fig. [3b that plots the energy distribution calculated 
using Eq. [4] at selected values kp along the arc of the un- 
reconstructed Fermi surface. To simulate experimental 
curves, the spectral intensity is modulated by /(ek — £f) 
at T — and broadened by an 'instrumental resolution' 
gaussian of width » |A. When symmetrized [8j], the 
energy distribution closely reproduces the apparent 
"pseudogap behavior" of P^S^CaC^Os+i [1]. 
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FIG. 3: a The approximately sinusoidal variation of the spec- 
tral weight in Fig. [2] (£ = 10 A) with angle 8, following the 
arc of the unreconstructed Fermi surface from Fig. [TJi. 6 is 
the angle between the line linking (ir, n) to (0, 0) and that 
linking (it, it) to the arc (inset), b Calculated energy distribu- 
tion intensity at kF along the unreconstructed Fermi surface 
at different 9. Red curves are calculations; black curves are 
"symmetrized" in the same manner as experimental data 

Quantum-oscillation experiments give a direct measure 
of ^4k,Q (Fig- lb) that is also subject to statistical vari- 
ations. The probability distribution is given by 



PA k = J J S(A k ^Q - A k )pQdQ x dQy. 



(5) 



Fig. fjji shows Eq. O for several £ ; the curves resemble 
the probability distributions caused by finite quasipar- 
ticle lifetime or other forms of disorder [l], [2(|, lead- 
ing to a qualitatively similar Landau-level broadening 
and quantum-oscillation damping. Figure 2)d shows pre- 
dicted damping factors for the fundamental quantum- 
oscillation frequency F ~ 650 T in YBa2Cu40s 0] and 
for the 2F harmonic. The correlation lengths £ « 10 and 
20 A that are best able to reproduce the k-dependent 
spectral weight of Cai.gNao.iCuC^C^ [7| in Fig.[5]cause 
too much damping to mimic the data of Ref. [1| ; a damp- 
ing factor ~ 10~ 3 is the typical cutoff for observable 
quantum oscillations, corresponding to £ w 90 A. An 
upper limit of £ < 400A comes from the absence of 
harmonics in the data, normally abundant in quasi-two 
dimensional metals i.e. 90 < £ < 400 A. The obser- 
vation of quantum oscillations in YBa2Cu306.5 [1| and 
YBa2Cu408 [3] therefore suggests a correlation length 
several times longer than that in ARPES experiments. 
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FIG. 4: a The >lk,q probability distribution for several val- 
ues of £. Grey lines mark the positions of the ideal 5-function 
Landau levels for a field of 80 T expected in this quasi-two di- 
mensional metal were long-range antiferromagnetic order es- 
tablished. These would furnish quantum oscillations periodic 
in 1/B, rich in harmonic content, b A plot of the correspond- 
ing damping factor versus £ for the fundamental frequency 
(black) and harmonic (red) resulting from the convolution of 
the probability distribution with the ^-function Landau levels. 

The differences in £ result from the contrasting con- 
ditions under which ARPES and quantum-oscillation 
experiments are performed. Whereas ARPES experi- 
ments are performed at tens of kelvins and zero magnetic 
field, quantum-oscillation experiments are made at liq- 
uid helium temperatures in the strongest fields available. 
Abundant neutron scattering data on Sr- and O-doped 
La2CuC>4 (as well as electron-doped cuprates) show an 
enhancement of antifcrromagnctism and/or its associ- 
ated correlation length at liquid 4 He temperatures un- 
der magnetic fields with values as large as £ ~ 
400 A reported 23 1. While similar neutron-scattering 
data are unavailable for the two systems (YBa2Cu3C>6.5 
and YBa2Cu40s) in which quantum oscillations are ob- 
served, alternative experimental techniques do provide 
evidence for field- induced antiferromagnetisrn [24| ■ These 
include nuclear magnetic resonance 25[ , muon spin rota- 
tion [26] and magnetic torque 27 1. 

In conclusion, we show that the seemingly contradic- 
tory results of ARPES and quantum oscillation experi- 
ments on the underdoped cuprates can be brought into 
mutual agreement by considering the slow antiferromag- 
netic fluctuations with a short correlation length detected 
by neutron-scattering experiments. Whereas the corre- 
lation length is very short (10 < £ < 20 A) in the nor- 
mal state at zero magnetic field, it becomes enhanced 
(90 ^ £ ^ 400 A) at liquid He temperatures in mag- 
netic fields due to a field-induced stabilization of anti- 
ferromagnetism [l6[. Short correlation lengths lead to a 
reduction in the spectral weight that is most pronounced 
close to (tt, 0) on the unreconstructed Fermi surface, giv- 
ing rise to the emergence of Fermi arcs near (tt, 7t) in a 
qualitatively similar fashion to a rf-wave order parame- 
ter. The effects described in this paper should also occur 
for other forms of ordering that lack long range period- 
icity including stripes and/or charge order Q. However, 
only antiferromagnetism is truly ubiquitous, giving rise 
to clear signatures in neutron-diffraction experiments in 



all underdoped cuprates investigated. 
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